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Abstract 

The double J/if>(DJ) production in direct photon-photon collision is investigated. 
It is found that the J /if) production rate in this process is of the same order of 
magnitude as those of previously discussed ones, which hints the dominant J/if> 
inclusive production process in the direct photon-photon collision is still not being 
touched. As DJ production process has a clear distinction from the single J/if> 
inclusive processes in the final states, experimental study on it would be helpful 
to clarify the validity of color-singlet model. 
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The quarkonium production and decays have long been taken as an ideal means in in- 
vestigating the nature of QCD and other new phenomena. Hence, to establish a proper 
theory which can precisely describe the heavy quarkonium production and decays is very 
necessary. The novel effective theory, the non-relativistic QCD(NRQCD) |J, is possibly the 
one to this aim, which is formulated from the first principles. To make a precise prediction 
for quarknoium production, however, at now, with only the NRQCD is not enough since the 
magnitude of the non-perturbative parameters in the theory are still unknown. Or in other 
words, results are always sensitive to the values of the input non-perturbative parameters, 
whereas the NRQCD can at most give the relative weights of these parameters in order of 
v 2 based on the "velocity scaling rules". The present situation in the Onium physics is that 
on one side the color-octet model @ stands as the most plausible approach, till now, in 
explaining the large transverse momentum ip{ip') production "anomaly"; on the other hand, 
the model meets some difficulties in confronting with other phenomena ||. Therefore, to 
what degree the color-octet mechanism plays the role in quarkonium production is still not 
clear and interesting. Either to find distinctive signatures of Color-octet or more strong 
arguments against it is currently an urgent task in this research realm. To this end lots of 
works have been done, however, unfortunately up to now there is still no definite answer. 

In recent years several new concepts on linear colliders aiming at providing collisions 
at the center-of-mass energy from hundreds GeV to multi-TeV with high luminosity are 
proposed and the feasibilities are pre-tested, such as JLC at KEK, TESLA at DESY and 
CLIC at CERN, etc. Theoretically, high-intensity photon beams may be obtained by the 
Compton back-scattering of laser light off the linac electron beams and realize the photon- 
photon collision with approximately the same luminosity as that of the e + e~ beams. Such 
a photon linear collider can have high energy up to the TeV order. Since the first discovery 
of J /ip in 1974, researches on quarknoium production at e + e~ colliders at various energies 
have been carried out in details |§-f|. However, studies on the photon-photon scattering are 
very limited and just begin [p|-[Tl[]. One of the difficulties one may meet in doing a complete 
calculation for two photon process beyond LO in a, the electromagnetic coupling constant, 
or a s , the strong coupling constant, is how to simplify the lengthy expressions in order to 
present them analytically. In this brief report we study the double quarkonium, the J/ip, 
production in direct photon-photon collision(of the resolved case see, e.g., ref. ||) at the 
order a 2 a 2 .. From physical point of view, in this process one may get rid of the inaccuracy 
induced by the uncertainties of the color-octet matrix elements and by the contributions 
from higher excited states as well, because they are much suppressed relative to the direct 
contribution in color- singlet. 

In 77 scattering, at leading order in a the J / ip is produced via the process 
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7 + 7^^/V + 7 



(1) 



However, since at the scale of heavy quark mass, the strong coupling constant is not too 
small, the process 

1 + 1 -*J/^+g (2) 



may compete with the pure electromagnetic process (p]) through the Color-Octet mechanism 
||. Here, we schematically use the J/i})^ to denote those evolved from the Color-Octet 
states. 

To go one order up in a s , one may still expect to get the same order of maginitude in J/ip 
production rate, because in this case the J / ip may be produced in color-singlet and therefore 
will get compensation for the a s suppression from the non-perturbative sector relative to the 
octet process @). What we proceed here, the DJ production in 77 collision, is a sub-category 
of the inclusive J/ip production process at order ot 2 a 2 s . That is 

7(*0 + 7(* 2 ) - J/i>(P) + J/^{P') , (3) 



as shown in Figure 1, which consists of twenty Feynman diagrams. Here, the momenta of 
the related particles are shown in the parentheses explicitly. Similar as the situation of B c 
production in photon-photon collision ||12|| , these twenty diagrams involved here can also be 
classified as two sub-groups, (A) and (B), and each of them are gauge invariant. 

The differential cross-section of process (|]) is obtained as 
da 262144a 2 a 2 7r|_R(0)| 4 9n 18/ , . A/ , 9 „ 9 , % , ,0 

Tt = -ty { rJ- W [mm ~ nl2m {t + u) + t u (t + u) ~ Sm t u (t + u) 

+ m 16 (1917t 2 + 3578tw + 1917m 2 ) - 8m 14 (242t 3 + 657t 2 w + 657tw 2 + 242m 3 ) 
+ 2m 12 (623t 4 + 2204t 3 w + 3298tV + 2204tw 3 + 623m 4 ) 

- 4m 10 (131t 5 + 563t 4 w + 1156tV + 1156*V + 563tw 4 + 131m 5 ) 

+ m 8 (143t 6 + 690t 5 w + 1899* V + 2640tV + 1899tV + 690tw 5 + 143u 6 ) 

- 4m 6 (6t 7 + 29t 6 u + 107tV + 212tV + 212tV + 107tV + 29tu 6 + 6u 7 ) 

+ 2m 4 (t 8 + At 7 u + 21tV + 68*V + 104tV + 68tV + 21tV + Atu 7 + u 8 )] . (4) 



In obtaining the above analytical expression, we start from the general Feynman rules and 
project the Charm- ant i- Charm pair into the S-wave vector Onium state in color- singlet. To 
manipulate the trace and matrix element square of those twenty diagrams, the computer 



algebra system MATHEMATICA is used with the help of the package FEYNCALC [13|. 
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The total cross section of the direct photon-photon collision can be obtained by convo- 
luting the differential cross-section with the photon distribution functions, like 

atotai = J dtdxidx 2 f 1 {xi)f 1 {x 2 ) < ^ L {xi,X2) , (5) 

where ^^-(xi, x 2 ) is the differential cross-section given in (f|); the f 7 (xi), i = 1,2, is the 
photon distribution with the fraction of the beam energy, 



1 — x H 4r(l — r) 

1 — x 



(6) 



Here, r = x/ (x m (l — x)), the normalization 



Xm ^rn ^ X m 2(1 -|- 



and 



X m = 5— COS - , (8) 

?7Jg 2 

where and £7 are the energies of electron beam and laser photon, respectively, and the 9 
is the angle between them. The energy fraction x of the photon is restricted in 

< x <^-, (9, 
and the maximum energy of the obtained photon beam depends on x m , i.e., 

Umax = : /" -Eft • (10) 

In our numerical calculation we take the optimum value of x m = 4.83 in avoiding the 
background e + e~ pair production from the laser and backscattered photon collision |T4 . 



Other values of input parameters used in our numerical calculation are: 
m c = 1.5 GeV, a = 1/128, a s (m c ) = 0.3, |R(0)| 2 = 0.8 GeV 3 , < >= O.OlGeV 3 . (11) 



With the above formulas and input parameters one can immediately get the total cross- 
sections. For instance, the J/ifj production cross-sections of the processes ([1]), (fj) and (|3|) 
are 2.54 fb, 1.92 fb and 2.15 fb, respectively with colliding energy being of 500 GeV. Since the 
projected linear colliders with luminosity of hundreds fb -1 per year, for example the TESLA, 
and the integrated total cross-sections increase with the colliding energy decreasing as shown 
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in Figure 2, we may have hundreds of events being observed in one year at colliding energy 
500 GeV or less. However, when the colliding energy goes up to ITeV the cross-sections for 
processes (1) - (3) would be three orders less than those of them at 500 GeV, respectively, 
hence, all of them would be unobservable. In Figure |2] the energy dependences of the total 
cross-sections are shown from about the energy of LEP II to the next generation of linear 
collider energy. It can be seen that the color-octet process is the smallest one all over this 
energy scope in the concerned three processes, though they are in the same order. 

Another thing we want to point out here is that the DJ production discussed in this work 
can not be simply explained by the quark fragmentation. With more details: we have the 
differential cross-section of 77 — > QQ as 



da —2Tca 2 eQ 



[6m 8 Q - tu{f 4~ it 2 ) — m^(3t 2 + Utu + 3u 2 



dt s 2 (rriQ — t) 2 (rriQ — u) 2 

+ m 2 Q (t 3 + 7t 2 u + 7tu 2 + u 3 )] , (12) 



and know the probability of Charm quark evolving to J/ip is of 1.6 x 10 4 via fragmentation 



| 15| . To a certain degree of approximation, one can get the DJ production rate through 
fragmentation mechanism by simply combining the Charm production rate with the universal 
fragmentation probability. Therefore, we have a cross-section of 1.1 x 10 _3 fb at energy 500 
GeV, which is much smaller than that of the complete calculation of the DJ process. The 
reason of this is that the fragmentation mechanism only applicable when the fragmenting 
quark possesses large momentum and small off-shellness. In the double Onium production 
this condition can not be fulfilled for both fragmenting quark at the same time. On the 
other hand, one can know from the above discussion that the cross-section for single J/ip 
production via fragmentation mechanism would be about one order larger than those of the 
processes (1) - (3) at 500 GeV. Furthermore, we get know that the process at order a 2 a 2 
for single J/ip inclusive production should be the dominant one in photon-photon collision, 
which seems to be overlooked. 

Considering there are remaining disagreements among three independent calculations on 
the leading order processes ([I]) and (^) P|-|TT|], we have done a check over them and find an 
agreement with refs. and ]TT[ on the expression for octet process, i.e., the eq.(7) of ref. 
0. And our replacement for obtaining the singlet process agrees with that in ref f9|. 

In conclusion, we have calculated the DJ production at photon colliders. It is found 
that at moderate energy of the next generation linear colliders there would be hundreds of 
events to be detected per year with the high projected luminosity. In addition, since the 
production rates of J/ip via the color-octet mechanism, electromagnetic process, and what 
discussed here are almost the same in the full scope of the colliding energy, to differentiate 
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the color-octet mechanism from the color-singlet one in these processes, experimentally one 
should detect not only the J/tp but also other final states. Relatively the processes (1) and 
(3) are more liable to be discerned than (2), because the signal of the latter is easily to be 
entangled with that of the single J/ip production process at the order a 2 a 2 s [16] . 



Like in the case of double n production at photon colliders, where the Pion wavefuctions 
and the photon-to-pion transition form factor may be detected [|T7|, the DJ process may pro- 
vide more information on the color-singlet description of quarkonium production considering 
that it gets less influences from the color-octet uncertainties and the higher excited states 
feeddown. 

In the end, with the expression (4) similar process in the soft bremsstrahlung photon 
collision can be easily discussed with the Weizacker- Williams approximation. 
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FIG. 1. The Feynman diagrams of the double J ftp production in 77 collision. 
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FIG. 2. The engergy dependence of the cross-sections. DJ: the double J/tp process; CS 
color-singlet process (H); CO: the color-octet process @. 
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